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Abstract 

A mechanistic, dynamic model was developed to predict infection of loquat fruit by conidia of Fusicladium eriobotryae, the 
causal agent of loquat scab. The model simulates scab infection periods and their severity through the sub-processes of 
spore dispersal, infection, and latency (i.e., the state variables); change from one state to the following one depends on 
environmental conditions and on processes described by mathematical equations. Equations were developed using 
published data on F. eriobotryae mycelium growth, conidial germination, infection, and conidial dispersion pattern. The 
model was then validated by comparing model output with three independent data sets. The model accurately predicts the 
occurrence and severity of infection periods as well as the progress of loquat scab incidence on fruit (with concordance 
correlation coefficients >0.95). Model output agreed with expert assessment of the disease severity in seven loquat- 
growing seasons. Use of the model for scheduling fungicide applications in loquat orchards may help optimise scab 
management and reduce fungicide applications. 
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Introduction 

Scab, caused by the plant-pathogenic fungus Fusicladium 
eriobotryae (Cavara) Sacc, is the main disease affecting loquat in 
Spain and in the whole Mediterranean basin [1,2]. The fungus 
affects young twigs, leaves and fruits, causing circular olive-colored 
spots that, on fruits, reduce their commercial value [1]. 
Fusicladium spp. are the anamorphic stages of the ascomycete 
genus Venturia but the sexual stage of F. eriobotryae has never 
been found in nature [2]. 

Although loquat scab is a well-known problem in the areas 
where loquat trees are cultivated, the biology of F. eriobotryae and 
the epidemiology of the disease have been seldom studied [1,3-8]. 
These studies have depicted F. eriobotryae as a highly rain- 
dependent pathogen that requires mild temperatures and long wet 
periods to infect loquat trees. 

Environmental requirements for infection and the dispersion 
patterns have been studied in detail for other Venturia spp., such 
as Venturia inaequalis [9-16], V. nashicola [17-19], V. pyrina 
[20-24], F. carpophilum [25,26], F. effusum [27-29], and F. 
oleagineum [30-35]. These studies have been used to elaborate 
epidemiological models for some of these pathogens including V. 
pyrina [36], V. nashicola[?<T\, V. inaequalis [38,39], F. oleagineum 
[40], and F. effusum [41]. For V. inaequalis, the use of 
epidemiological models to schedule fungicide applications has 



reduced the number of treatments [42-45]. To date, no 
epidemiological model has been developed for F. eriobotryae . 

Disease modelling is an important step towards the implemen- 
tation of sustainable agriculture [46,47]. Since the 1990 s, modern 
crop production has focused on the implementation of less 
intensive systems with reduced inputs of fertilizers and pesticides, 
and reduced use of natural resources [46] . Sustainable agriculture 
has its roots in Integrated Pest Management (IPM) [48]. IPM 
concepts originated as a reaction to the disruption of agro- 
ecosystems caused by massive applications of broad-spectrum 
pesticides in the middle of the last century [46] and also because of 
concern about the effects of excessive pesticide use on human 
health [49]. 

In Europe, the implementation of IPM has been legislatively 
mandated in recent years because of Directive 2009/128/CE 
regarding sustainable use of pesticides. Among other actions, the 
Directive encourages EU Member States to promote low 
pesticide-input pest control and the implementation of tools for 
pest monitoring and decision making, as well as advisory services 
(Art. 14 of the Directive). De facto, the "sustainable use" directive 
has made IPM mandatory in European agriculture as of 20 1 4. As 
a consequence, there is an increased interest in the development 
and use of plant disease models to improve the timing of pesticide 
applications and to thus limit unnecessary treatments [46,50,51]. 
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Our aims in this paper were (i) to develop a mechanistic, 
dynamic model to predict infection of loquat fruit by the scab 
fungus F. eriobotryae, and (ii) to evaluate the model against three 
independent data sets. The model was elaborated based on the 
principles of "systems analysis" [52,53] and by using recent data 
on the biology and epidemiology of F. eriobotryae obtained under 
environmentally controlled and field conditions [1,6,7]. 

Model Development 

Based on the available information [1,3—8], the life cycle off. 
eriobotryae under the Mediterranean climate is described in 
Figure 1. The fungus oversummers in lesions on branches and 
leaves and on mummified fruits that remain in the tree after 
harvest; during summer, high temperatures and low humidity may 
prevent sporulation on these lesions. Under favorable conditions in 
the fall, the conidia produced by the oversummering lesions serve 
as the primary inoculum and infect young leaves or loquat fruits. 
Conidia are dispersed by splashing rain to nearby fruits and leaves; 
with suitable temperature and wetness, conidia germinate and 
penetrate the tissue, probably direcdy through the cuticle or 
through stomata. Once infection has occurred and if the 
temperature is favorable, the fungus grows under the cuticle; 
conidiophores then erupt through the cuticle and produce new 
conidia. These conidia cause secondary infections during the 
entire fruiting season as long as rains disperse them and as long as 
temperature and wetness duration permit conidial germination, 
infection, and lesion growth. 

Model description 

The relational diagram of the model for loquat fruit infection by 
F. eriobotryae is shown in Figure 2, and the acronyms are 
explained in Table 1 . The time step of the model is 1 hour. 



The model starts at fruit set and ends at harvest because fruits 
are assumed to be always susceptible to infection. The model 
considers the lesions from the previous season on branches, old 
leaves, and mummified fruits as the sources of primary inoculum. 
Because the abundance of these lesions in an orchard may vary 
depending on several conditions-on, for instance, the level of 
disease or the fungicide treatments in the previous season-and 
because it is difficult to quantify these lesions, the model assumes 
that oversummered forms are present in the orchard and that they 
hold conidia at fruit set and onwards. 

The model considers that any measurable rain (i.e., R^0.2 mm 
in 1 hour) causes dispersal and deposition of conidia on loquat fruit 
[7] and triggers an infection process that potentially ends with the 
appearance of scab symptoms. Each site on the fruit that is 
occupied by a conidium or conidia is considered a potential 
infection site and is referred to as a lesion unit (LU). During the 
infection process, infection on any LU can fail because conidia 
may fail to germinate or may germinate but then die because of 
unfavorable conditions. Therefore, the proportion of LUs that 
become scabbed at the end of the infection process may be less 
than that occupied by splashing conidia at the beginning of the 
process. 

The model predicts the progress of infection on single LUs, 
which are the surface unit of the fruit which can become occupied 
by a scab lesion. This approach is related to the concept of 
"carrying capacity". In ecology, the carrying capacity is 
interpreted broadly as the maximum population size that any 
area of land or water can sustain [54,55] . In plant pathology, the 
host's carrying capacity for disease is the maximum possible 
number of lesions that a plant (or an organ) can hold [56] . The 
carrying capacity is a common concept in plant disease modeling 
[57-60]. In the model described here, a LU is initially healthy 
(LUH) but then becomes occupied by: ungerminated conidia 




Figure 1. Disease cycle of loquat scab caused by Fusicladium eriobotryae. 

doi:1 0.1 371 /journal.pone.01 07547.g001 
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Table 1. List of variables used in the model. 



Acronym 


Description 


Unit 


T 


Air temperature 




DU 
tin 


Relative humidity 


70 


Q 

n 


Rainfall 


mm 


VrU 


Vapour pressure deficit 


hp- 

urd. 


\A/n 

VvU 


Wetness duration 


hours 


Tea 


Temperature equivalent 


c 


Lun 


Unit of loquat fruit surface without conidia of F. criobotryoc 


iNumDer (u— \ ) 


LUUC 


Unit of loquat fruit surface with ungerminated conidia of F. criobotryoc 




LUGC 


Unit of loquat fruit surface with germinated conidia of F. criobotryoc 


INumDer — i y 


LULI 


Unit of loquat fruit surface with latent infection by F. criobotryoc 


M i imkor /fl 1 ^ 

iNumuer \*J~ i ) 


LUVI 


Unit of loquat fruit surface with visible scab lesions 


iNumuer \ ) 


SZPD 


Cumulated conidial germination 


Kill mknr i(\ 1 \ 

iNumuer \ ) 


INF 


Cumulated infection 


Number (0-1) 


SUR 


Cumulated conidial survival 


Number (0-1) 


GER' 


Germination rate (first derivative of GERM) 


Number (0-1) 


INF' 


Infection rate (first derivative of INF) 


Number (0-1) 


SUR' 


Survival rate (first derivative of SUR) 


Number (0-1) 


C 


Correction factor 


Number (0-1) 


DD 


Degree days 


Number 



doi:1 0.1 371 /Journal.pone.01 07547.t001 




Figure 2. Relational diagram showing how the model simulates infection by Fusicladium eriobotryae. Legend: boxes are state variables; 
line arrows show fluxes and direction of changes from a state variable to the next one; valves define rates regulating these fluxes; diamonds show 
switches (i.e., conditions that open or close a flux); circles crossed by a line show parameters and external variables; dotted arrows show fluxes and 
direction of information from external variables or parameters to rates or intermediate variables; circles are intermediate variables. See Table 1 for 
acronym explanation. 
doi:1 0.1 371 /journal.pone.01 07547.g002 
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(LUUC) at the time of conidial dispersal; germinated conidia 
(LUGC) at the time of conidial germination; latent infection after 
penetration (i.e., hyphae are invading the fruit cuticle; LULI); and 
visible and sporulating scab lesions at the end of latency (LUVI). 
Both LUUC and LUGC can fail to progress if ungerminated or 
germinated conidia die; these LUs then return to being LUHs 
because they can start a new infection process whenever new 
conidia are splashed on them. 

At any dispersal event on hour h, the model considers that 
LUUCh = 1 . The rate at which LUUC^ advances to LUGCj, 
depends on a germination rate (GER'), and the rate at which 
LUGC h advances to LULI h depends on an infection rate (INF') 
(Figure 2). Both GER' and INF' are influenced by temperature (T 
in °C) and wetness duration (WD, in hours) (i.e., free water on the 
surface of the loquat fruit) caused by either rain or dew. Fruit 
surfaces are assumed to be wet on any hour when Rh>0 mm, or 
RH h >89%, or VPD h <l, where VPD is the vapour pressure 
deficit (in hPa) calculated using T h and RH h , following Buck [61]. 
The rate at which LUUC^ and LUGCj, returns to LUH^ depends 
on a survival rate (SUR'), which depends in turn on the length of 
the dry period (DP), i.e., the number of hours with no wetness on 
the fruit surface (Figure 2). 

GER', INF', and SUR' are calculated at hourly intervals by 
using the first derivative of the equations described in Gonzalez- 
Dominguez et al. [6] in the form: 

GER' = 

116.249 x Teq 4347 x (1 - Teq 2M2 ) x (!) 
3.215x e <- 0 - 376x ^)x e [- 8 ' 551x ^°' 376xra) ] 

INF' = 

4.961 x Teq ljm x (1 - Teq) 0Jn x (2) 
0.409 x e - 008 ™x e L J 



L UGC h +LUUC h +L ULI h + LUH h = \ 

where C is a correction factor (C = 1 — LULI/,) 

Any infection period triggered by a conidial dispersal event ends 
when no viable conidia are present on any LUs, exactiy when 
LUUC^O.Ol. An example of model output for a single infection 
period is shown in Figure 3. 

The model considers that any further rain event causes further 
dispersal and deposition of conidia if >5 hours have passed after 
the previous dispersal event. This is the time required by a lesion 
to produce new conidia. 

Model output 

The model output consists of: (i) the available inoculum on fruits 
(i.e., the frequency of LUs with ungerminated conidia on each day) 
as a measure of the potential for infection to occur; (ii) the 
dynamics of LULI for each infection process; and (iii) the seasonal 
dynamics of the accumulated values of LULI (ZLULI) as an 
estimate of the disease in the orchard. 

Examples of model output for the 2011 and 2012 loquat 
growing seasons are shown in Figures 4 and 5, respectively. The 
output is based on the weather data registered by a weather station 
of the Regional Agrometerological Service (http://riegos.ivia.es/) 
located in Callosa d'En Sarria, Alicante Province, southeastern 
Spain. 

Model Validation 

Three data sets were used to validate the model: (i) incidence of 
affected fruits in a loquat orchard during growing seasons 20 1 1 
and 20 1 2; (ii) disease occurrence on loquat fruits in single-exposure 
experiments in 2013; and (iii) expert assessment of the disease 
severity in seven loquat growing seasons. 

To operate the model, hourly values of air temperature (T, °C), 
relative humidity (RH, %), and total rainfall (R, mm) were 
registered by the weather station of Callosa d'En Sarria, which is 
<3.5 km from the orchards considered for validation. 



SUR' 



0.165 
DP 



(3) 



where: Teq is the temperature equivalent in the form 
Teq=(T—T m i n )/(T max — T m in) where: T is the temperature 
regime, T mm = 0°C and T max — 35°C in equation (1), and 
T min = 0°C and T max — 25°C in equation (2); WD = number of 
consecutive hours with wetness; DP = number of consecutive 
hours with no wetness. When DP = 0, SUR' = 1 
At any time of the infection progress (i): 



LUGC h = GER '. x 0 - J2 SC//? ' ! ') x C ' 



LUUC h = J2 (i-LUGC h )i x (1 - SUR'i) x C, 



LULI h = ^2 (INF'i 



Predicted vs. observed disease incidence in orchards 

In data set (i), observations were carried out in a loquat orchard 
in Callosa d'En Sarria, Alicante Province, southeastern Spain. 
Details on these data have been previously published [7]. Briefly, 
fruits from four shoots of each of 46 loquat trees were assessed 
weekly, and disease incidence was expressed as the percentage of 
fruits with scab symptoms. The disease incidence was lower in 
201 1 than 2012, with 27.3% and 97.6% of fruits affected by loquat 
scab at harvest, respectively [7]. This difference in disease 
incidence may be related to the fact that the orchard was treated 
with fungicides for scab control in 2010 but not in 201 1 or 2012. 
Given that the inoculum sources for fruit infection in 2011 was 
very low because of effective disease control in 2010, a correction 
factor for LUUC was applied for the infection processes initiated 
in January 201 1, i.e., LUUC = 0.1 instead of = 1 in January 201 1. 

Model validation was performed by comparing ZLULI with 
observed data of disease incidence. Because there is a time lag (i.e., 
a latency period) between the predicted disease (as 2LULI) and 
the disease incidence estimated in the orchard (DI), DI was shifted 
back by one latency period for comparison between predicted and 
observed disease. Sanchez-Torres et al. [1] observed a latency 
period of 2 1 days at a constant temperature of 20°C, which is a 
degree-day accumulation (DD base 0°C) of 420. Therefore, DI 
was shifted back by either 21 days or 420 DD. To calculate the 
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Figure 3. Dynamics of lesion units (LUs) during an infection period of Fusidadium eriobotryae. The graph shows the relative frequency of 
LUs occupied by ungerminated conidia (LUUC, in green), germinated conidia (LUGC, in red), and latent infections (LULI in blue). Blue bars at the top 
indicate hours with free water on the fruit surface. An infection period starts when a rain event splashes conidia on LUs and ends when no viable 
conidia are present on any LUs, i.e., when LUUCsO.01. 
doi:1 0.1 371 /journal.pone.01 07547.g003 



DD, the average temperature of each day was considered with 
base temperature of 0°C. 

Predicted vs. observed disease incidence in single- 
exposure experiments 

In data set (ii), data were collected in an abandoned loquat 
orchard in Callosa d'En Sarria from 4 February to 15 April 2013. 
On 25 January, 200 random shoots bearing fruits were covered 
with water-resistant paper bags (one shoot per bag) to prevent 
deposition of rain-splashed conidia. On 4 February, 10 random 
bags were opened to receive splashed inoculum; after seven 
additional days, the bags were closed again. Ten other randomly 
selected bags were opened on 1 1 February and closed again 7 days 
later. This operation was repeated until nine groups of shoots had 
been sequentially exposed to rain. At the end of the experiment (15 
April 2013), disease incidence (percentage of fruits affected by 
loquat scab) and severity were assessed in each group of shoots. 
Disease severity refers to the percentage of fruit area covered by 
scab lesions and was measured as described by Gonzalez- 
Dominguez et al. [8]. 

Model validation was performed by comparing the model 
output in the week when a group of shoots was exposed to 
splashing rain with final disease severity in that group. 

Expert assessment 

For data set (iii), Esteve Soler (technical advisor of the 
'Cooperativa Agricola de Callosa d'En Sarria') was asked to 
provide a subjective estimate of the severity (low, medium, or high) 
of loquat scab in the area for eight growing seasons (from 2005/ 
2006 to 2012/2013). Mr. Soler's estimates were based on his 



extensive experience in managing loquat orchards, on his scouting 
activities in the orchards of the cooperative, and on the number of 
fungicide treatments that were required to control the disease in 
the area. 

For each season, the model was operated from 1 November to 
3 1 March, and the numbers of disease outbreaks predicted by the 
model were counted. A disease outbreak was defined as ELULI> 
0. 1 in 1 day, when no outbreaks were predicted in the previous 5 
days. Average and standard error of the number of predicted 
outbreaks were calculated for each category (low, medium, or 
high) of scab severity derived from the expert assessment. 

Data analysis 

Linear regression was used to compare the predicted and 
observed data of data sets (i) and (ii). To make data homogeneous, 
ELULI values at the time of each disease assessment in the 
orchards were rescaled to the ELULI at the end of the season; 
disease incidence was also rescaled to the final disease incidence. A 
i-test was used to test the null hypotheses that "a" (intercept of 
regression line) was equal to 0 and that "b" (slope of regression 
line) was equal to 1 [62]. The distribution of residuals of predicted 
versus observed values was examined to evaluate the goodness-of- 
fit. The concordance correlation coefficient (CCC) was calculated 
as a measure of model accuracy [63]; CCC is the product of two 
terms: the Pearson product-moment correlation coefficient 
between observed and predicted values and the coefficient Cb 
(bias estimation factor), which is an indication of the difference 
between the best fitting line and the perfect agreement line 
(CCC = 1 indicates perfect agreement). The following indexes of 
goodness-of-fit were also calculated [64]: NS model-efficacy 
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Leaf wetness (h) tmm Rain (mm) Temperature (°C) 




Figure 4. Weather data and model output in 201 1. A: daily weather data; B: predicted frequency (%) of lesion units (LUs) with ungerminated 
conidia; C: predicted increase of LUs with latent infections (LULls) for each infection period (arrows represent clusters of infection periods, clustering is 
based on an interval of at least 5 days between the beginning of two consecutive clusters); D: predicted seasonal dynamics of the cumulative values 
of LULI (ELULI). 

doi:1 0.1 371 /journal.pone.01 07547.g004 



coefficient, which is the ratio of the mean square error to the 
variance in the observed data, subtracted from unity (when the 
error is zero, NS = 1, and the provides a perfect fit); the W index of 
agreement which is the ratio between mean square error and total 
potential error (W = 1 represents a perfect fit); model efficiency 
(EF) which is a dimensionless coefficient that takes into account 
both the index of disagreement and the variance of the observed 
values (when EF increases toward 1, the fit increases); and the 
coefficient of residual mass (CRM) which is a measure of the 
tendency of the to overestimate or underestimate the observed 
values (a negative CRM indicates a tendency of the model toward 
overestimation). 

For data set (iii), a one-way analysis of variance (ANOVA) was 
performed to determine whether the numbers of outbreaks 
predicted by the model in each category of loquat scab severity 
defined by the expert (i.e., low, medium, or high) were significantly 
different from one another. 



Results of Model Validation 

Predicted vs. observed disease incidence in orchards 

In 2011 between 1 January (fruit set) and 23 May (harvest), 
257.6 mm of rain fell, distributed in three main periods: the last 
week of January, the second week of March (with 64.8 mm of rain 
in 1 day), and the last 2 weeks of April (with daily temperature > 
15°C) (Fig. 4A). According to the model, a total of 33 infection 
periods were triggered by these rain events, and the first was on 1 7 
January (Fig. 4B and 4C). In the analysis of this model output, 
infection periods were clustered in "infection clusters" based on an 
interval of a minimum of 5 days elapsed between the beginning of 
two consecutive infection clusters (i.e., the protection provided by 
a copper-based fungicide application as described in [65]); 
therefore, there were 10 infection clusters in the considered 
period. SLULI began to increase from mid-January to mid- 
February (with three infection clusters), but three infection clusters 
in March resulted in a substantial increase in ZLULI to 0.5; 
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Rain (mm) 



-Temperature (°C) 




Figure 5. Weather data and model output in 2012. A: daily weather data; B: predicted frequency (%) of lesion units (LUs) with ungerminated 
conidia; C: predicted increase of LUs with latent infections (LULIs) for each infection period (arrows represent clusters of infection periods, clustering is 
based on an interval of at least 5 days between the beginning of two consecutive clusters); D: predicted seasonal dynamics of the cumulative values 
of LULI (ELULI). 

doi:1 0.1 371 /journal.pone.01 07547.g005 



March had 12 infections periods, and the repeated and abundant 
rain events provided >18 h of wetness on most days (Fig. 4). From 
mid- April to the end of the considered period, a constant increase 
in ZLULI was associated with abundant rain events and 
increasing temperature, which triggered four infection clusters 
(Fig. 4). 

In 2012, although the total volume of rain that fell from 15 
December to 10 May was similar (255.8 mm) to that in 2011, 
there were fewer rain events. The model predicted 20 infection 
periods that were grouped into 12 infection clusters (Fig. 5B and 
5C). In 2012, rainy periods were separated by dry periods; from 
the end of January to mid-March, dry periods caused no 
substantial infection to develop (Fig. 5C). Therefore, there were 
two main periods of ELJJLl increase: the last half of January and 
from the end of March to May (Fig. 5D). 

Goodness-of-fit of predicted (ELULI) versus observed data 
(loquat scab incidence) was greater when a fixed period of 2 1 days 
was considered for the latency. In this case, values of R 2 , CCC, r, 



Cb, NS, W, and EF were >0.95 (Table 2). However, when a 
latency period of 420 DD was considered the values of R 2 and 
CCC were <0.88, and values of model efficacy (NS) and model 
efficiency (EF) were 0.75 (Table 2) as a consequence of the high 
dispersion of residues in 2012 (Figure 6). The model slightly 
overestimated scab incidence when a latency of 2 1 days was used 
(CRM = —0.009) and underestimated scab incidence when DD 
were used (CRM = 0.182) (Table 2; Figure 6). For both latency 
options, the regression equations of predicted versus observed data 
had slopes and intercepts that were not significantly different from 
1 and 0, respectively. 

Predicted vs. observed disease incidence in single- 
exposure experiments 

From 4 February to 15 April 2013, the model predicted 15 
loquat scab infection periods but disease outbreaks were substan- 
tial (i.e., they resulted in a >10% increase in severity) in only two 
exposure periods. In these two cases, LULI values were >0.1; 
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when there were no or light outbreaks, LULI values were <0.06 
(Figure 7). The goodness-of-fit of predicted versus observed for 
data set (ii) (Table 2) provided values >0.97 for R 2 , CCC, r, Cb, 
NS, W, and EF. Although the slope was not significantly different 
from 1, the intercept was different from 0 at P — 0.02 (Table 2). 
The negative value of CRM indicated that the model somewhat 
overestimated disease, mainly when observed disease severity was 
low (Figure 7). 
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Expert assessment 

The loquat scab epidemics that occurred in the eight seasons of 
data set (iii) were considered by the expert to be of low, medium, 
or high severity in two, three, and three seasons, respectively. The 
number of outbreaks predicted by the model ranged from 4 to 1 7 
among the eight seasons; in average, 8.5±0.5 outbreaks were 
predicted for years with low value of loquat scab severity, 1 0 ± 3 for 
year with medium value and 1 2 ±2.9 for years with high value. 
Although the average number of outbreaks predicted by the model 
increased as the expert assessment of disease severity increased, the 
number of predicted epidemics did not significantly differ among 
the severity categories (P = 0.71). 
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Discussion 

In this work, a dynamic model was developed to predict 
infection of loquat fruits by conidia of F. eriobolryae. The model 
uses a mechanistic approach to describe the infection process 
[53,66,67]: the model splits the disease cycle of F. eriobotryae into 
different state variables, which change from one state to the 
following state based on rate variables or switches that depend on 
environmental conditions by means of mathematical equations. 
The mathematical equations were developed using published data 
on F. eriobotryae conidial dispersion patterns [7] and on F. 
eriobotryae growth, conidial germination, and infection under 
different environmental conditions [1,6]. In the absence of precise 
information, assumptions were made based on available knowl- 
edge. 

Model validation showed that the model correctly predicted the 
occurrence of infection periods and the severity of any infection 
period, as demonstrated by the goodness-of-fit for the data 
collected on fruits exposed to single rainy periods. Because the 
purpose of the model is to be part of a warning system for loquat 
scab management, the ability to correctly predict infection periods 
is crucial. Accuracy of the model was also confirmed by the 
comparison of model output with expert assessment. Even though 
the numbers of predicted outbreaks did not differ among seasons 
that the expert had categorized as having low, medium, or high 
disease severity, the number of predicted outbreaks increased with 
increases in assessed disease severity. 

For model validation, the latency period required for the 
appearance of scab was expressed as a fixed number of days or of 
degree-days (DD) based on results from Sanchez-Torres et al. [1]. 
Goodness-of-fit of model prediction was overall better using a fixed 
period of 21 days instead of 420 DD. In particular, the model 
underestimated the disease in the early season of 2012 when DD 
were used. The underestimation was probably caused by low 
temperatures in that period, which delayed DD accumulation. 
This result is questionable, because the physiological development 
of fungi is usually more closely related to DD than to calendar days 
[68,69]. In this work, DD was fixed based on the latency period 
observed in loquat plants kept at the optimal temperature for F. 
eriobotryae development, i.e., 21 days at 20°C [1]. Therefore, the 
DD value used in this study did not account for the non-linear 
response oi F. eriobotryae growth to temperatures between 5 and 
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Figure 6. Comparison between model output and scab observed on loquat fruit in southeastern Spain. (A) data from 201 1 and (B) data 
from 2012. Blue lines represent the rescaled infection predicted by the model as the seasonal summation of the lesion units with latent infections 
(ELULI). Points represent rescaled incidence of loquat fruit with scab observed in the orchards; rescaled incidence is shifted back by 21 days (red 
points) or 420 DD (base O 'C, green points) to account for the latency period, i.e., the time elapsed between infection and visible symptoms in the 
form of sporulating scab lesions. 
doi:1 0.1 371 /journal.pone.01 07547.g006 



30°C [6]. If a function for predicting the appearance of scab 
symptoms is needed in the model, such a function should be 
temperature dependent, as it is hi models for V. nashicola [37] and 
F. oleagineum [40]. Salerno et al. [4] repeatedly exposed potted 
loquat plants under the canopies of affected trees for 3 days and 
then incubated these plants under a roof until the appearance of 
symptoms. Scab appeared in 11 to 26 days at temperatures 
ranging from 1 1.4 to 17°C (with a DD range of 157 to 340) and 
after >220 days at temperatures >20°C. Ptskialadze [5] found 
scab symptoms on both leaves and fruits 34 and 16 days after 
infection at 1-4°C and 21-25°C, respectively. Even though the 
calculation of latency can be improved, the model error in 
predicting disease onset due to a fix latency period may not reduce 
the ability of the model to correctly predict infection periods or 
reduce the value of the model for timing fungicides applications. 

The model capitalized on recent research concerning loquat 
scab [1,6,7]. These studies have considered most of the 
components of the disease cycle, including dispersion of conidia, 
infection, incubation, and latency. Nevertheless, other components 
should be elucidated to improve our knowledge and thus to 
improve the model [66]. Currently, the model assumes that 
inoculum sources are always present in scab-affected loquat 
orchards and that viable F. eriobotryae inoculum is always present 



at fruit set (i.e., when the model begins operating) and beyond. 
Salerno et al. [4] found that lesions appear in autumn on leaves 
that were infected the previous spring, and Prota [3] found that 
the lesions appearing in autumn produce conidia for 5 to 6 months 
and that those viable conidia are present all year long. These 
observations were carried out in Sicily and Sardinia, respectively 
(i.e., under a Mediterranean climate); therefore, the model 
assumptions seem plausible. The assumptions that inoculum 
sources and viable conidia are always present in scab-affected 
loquat orchards are both precautionary because they can lead to 
over prediction of infection (which would occur if weather 
conditions were suitable for infection but no viable conidia were 
available) and thus to unnecessary applications of fungicides or 
other disease management measures. Because unnecessary fungi- 
cide applications entail costs for growers, consumers, and the 
environment [51], the model should be expanded to include the 
oversummering and availability of conidia. 

With respect to oversummering, modeling the dormant stage of 
fungal pathogens is challenging [66], and the dormant stage has 
therefore been included in only a few models [70-74]. For this 
purpose, two key aspects must be addressed: (i) the inoculum dose 
(i.e., the quantity of inoculum that oversummers), which depends 
on the severity of the disease in each orchard at the end of the 
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Figure 7. Comparison between model output and scab on loquat fruit in single-exposure experiments. Experiments were carried out in 
a loquat orchard in southeastern Spain in 2013. Observed data (X axis) are expressed as the rescaled disease severity in 1 1 groups of fruits that were 
exposed (for 7-day-long moving periods) to splashing rain in a severely affected orchard; model output (Y axis) is expressed as the summation of the 
lesion units with latent infections (ELULI) in the exposure period. 
doi:1 0.1 371 /journal.pone.01 07547.g007 



previous season; and (ii) the time when the primary inoculum 
begins to be available for infection. In other models, the inoculum 
dose was directly measured in the field [70,74] or broadly 
estimated as low/high disease pressure [72]. In our case, 
incorporation into the model of the specific farmer's assessment 
of the disease severity in the previous season may represent useful 
information regarding the potential primary inoculum dose. 
Modeling the sporulation patterns of F. eriobolryae may make it 
possible to estimate the available inoculum at each infection 
period. This estimation may consequently improve the ability of 
the model to predict the severity of each infection period. To 
account for the presence of inoculum in a model for V. inaequalis, 
Xu et al. [39] assumed a minimum interval of 7 hours between 
two successive infection processes to allow lesions to recover and 
sporulate, even though this approximation could introduce errors, 
because sporulation is highly dependent on temperature and RH 
[14]. 

Even without the above possible improvements, the present 
model can contribute to the practical control of loquat scab. The 
underutilization of disease predictive systems by farmers have been 
broadly discussed [46,47,51,75,76]. Rossi et al. [53] summarized 
the steps necessary for the practical implementation of a model as: 
(i) develop a computerized version of the model; (ii) create a 
network of agro-meteorological stations for collecting weather 
data; (iii) design a strategy for decision-making based on the model 
output; (iv) develop tools for supporting decision-making (e.g., 
decision support systems or disease warning systems); and (v) build 
user's confidence in the model by demonstrating the advantages of 
its use in comparison with the current options. Efforts devoted to 
the last three steps are crucial for the future applicability of the 
model [53] and requires a deep knowledge of the cultural context 
in which the model will be delivered, the farmers' perception of 
risk, and the current management of the disease [47]. 

In the main loquat cultivation areas of Spain, the regional plant 
protection services use the Mills-Laplante tables [77], which were 



developed to control apple scab, to estimate the risk of infection by 
F. eriobotryae [78]. Researchers have indicated that the Mills- 
Laplante tables over-predict the number of infections for apple 
scab [37,79]. That the tables could over-predict the number of 
loquat scab infections has also been discussed, because the conidia 
of F. eriobotryae require longer times for leaf infection than those 
described by the Mills-Laplante tables for V. inaequalis and 
because the temperature range in which F. eriobolryae infection 
occurs is quite different [6,7]. Thus, the present model represents 
an improvement in loquat scab management, i.e., it should 
optimise scab management by helping loquat growers to schedule 
and probably to reduce fungicide applications. 

The long-term existence of a warning system for loquat scab 
monitoring in Spain [80] may facilitate the implementation of the 
model developed in this area because i) extension agents and 
advisors are familiar with the use and interpretation of epidemi- 
ological models, and ii) loquat farmers are accustomed to 
considering the concept of "infection risk" when scheduling 
fungicide applications. 

Because model building is "a never-ending story" [47,81], 
researchers will likely continue to improve the loquat scab model 
described here. As discussed in this manuscript, it will be necessary 
to define a relationship between model output and infection 
severity so as to identify appropiate thresholds for deciding when 
the treatments are needed [53,63]. 
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